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ABSTRACT 



We developed a new method for full kinematical reconstruction of the tt system near its 
threshold at future linear e + e~ colliders. In the core of the method lies likelihood fitting which is 
designed to improve measurement accuracies of the kinematical variables that specify the final 
states resulting from tt decays. The improvement is demonstrated by applying this method to 
a Monte-Carlo tt sample generated with various experimental effects including beamstrahlung, 
finite acceptance and resolution of the detector system, etc. A possible application of this 
method and its expected impact are also discussed. 
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1 Introduction 



The discovery of the top quark JTJj at Tevatron has completed the standard-model (SM) list 
of matter fermions. In spite of its subsequent studies thereat, however, our knowledge on 
its properties is still far below the level we reached for the lighter matter fermions. A next- 
generation e + e~ linear collider such as JLC []2[, having a facet as a top-quark factory, is expected 
to allow us to measure top quark's properties with unprecedented precision, thereby improving 
this situation dramatically. Such precision measurements may shed light on the electroweak- 
symmetry-breaking mechanism or hint beyond-the-SM physics or both. 

Being aware of the opportunities provided by the linear collider, a number of authors 
have so far taken up this subject and performed interesting analyses 0, |], ^, ||. These anal- 
yses can be classified into two categories, i.e., those near the ti threshold, mainly focused on 
physics contained in the threshold enhancement factor, and those in open-top region, searching 
for anomalies in production and decay vertices, both of which play important roles and com- 
plement each other. In order to thoroughly carry out such analyses for real data, we need a 
sophisticated method to kinematically reconstruct events as efficiently and as precisely as pos- 
sible. This is, however, highly non-trivial in practice, due to finite detector resolutions, possible 
missing neutrinos in the final states, and various background contributions. We thus need to 
further explore the potential of the e + e~ linear collider and extend the past studies 0, §, |J to 
make maximum use of the linear collider's advantages: clean experimental environment, well- 
defined initial state, availability of highly polarized electron beam, possibility of full parton-level 
reconstruction of final states, etc. 

As mentioned above, feasibility studies on form factor measurements are mainly done in 
the open-top region g |S], |TDfl . It has, however, been pointed out that form factor measurements 
in the ti threshold region have some favorable features such as availability of well-controlled 
highly polarized top sample, no need for transformation to t or i rest frames because both of 



the t and i are nearly at rest JTTJ, |I2 . Besides, as far as the decay form factor measurements 
of an on-shell top quark are concerned, the center-of-mass energy does not matter. Motivated 
by these observations, we decided to study feasibility of form factor measurements at the ti 
threshold. From an experimental point of view, the top quark physics at the e + e~ linear 
collider is expected to commence near the ti threshold and analysis techniques developed there 
are expected to be easily extended to the open-top region. 

In this paper, we thus aim at state-of-the-art analysis method such as a new method to 
reconstruct kinematically the ti system near its threshold in e + e~ annihilation. The content will 
be divided into six parts. In Sec. 2 we briefly review our simulation framework. Sec. 3 is devoted 
to top quark reconstruction in the lepton-plus-4-jet mode, where two subsections recapitulate 
basic strategy and procedure, respectively. In Sec. 4 we explain our kinematical reconstruction 
method which is designed to improve measurement accuracies of the kinematical variables that 
specify the final states resulting from ti decays. Then we discuss a possible application of this 
method and its expected impact in Sec. 5. Finally, Sec. 6 summarizes our results and concludes 
this paper. 



2 Framework of Analysis 



For Monte-Carlo-simulation studies of tt productions and decays, we developed an event gen- 
erator that is now included in physsim-2001a ||17||, where the amplitude calculation and phase 



space integration are performed with HELAS/BASES [|T8| , |19[ and parton 4-momenta of an event 
are generated by SPRING |19|| . In the amplitude calculation, initial state radiation (ISR) as well 
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as S- and P-wave QCD corrections to the ti system ]T3, |1(J are taken into account. Parton 
showering and hadronization are carried out using JETSET 7 . 4 ||20| with final-state tau leptons 
treated by TAUOLA [[?!]] in order to handle their polarizations properly. 

In this study, the top-quark (pole) mass is assumed to be 175 GeV and the nominal center- 
of-mass energy is set at 2 GeV-above the IS* resonance of the ti bound states. This energy is 
known to be suitable for measurements of various properties of the ti system at threshold |7j]. 
We will assume an electron-beam polarization of 80% in what follows. Effects of natural beam- 
energy spread and beamstrahlung are taken into account according to the prescription given in 
|?J , where the details of the beam parameters are also described. We have assumed no crossing 
angle between the electron and the positron beams and ignored the transverse component of 
the initial state radiation. Consequently, the ti system in our Monte-Carlo sample has no 
transverse momentum. Under these conditions we expect 40k ti events for \00fb~ 1 . 

The generated Monte-Carlo ti events were passed to a detector simulator (JSF Quick 
Simulator |22|]) which incorporates the ACFA-JLC study parameters (see Table. [I]). The 
quick simulator created vertex-detector hits, smeared charged-track parameters in the central 
tracker with parameter correlation properly taken into account, and simulated calorimeter 
signals as from individual segments, thereby allowing realistic simulation of cluster overlapping. 
It should also be noted that track-cluster matching was performed to achieve the best energy- 
flow measurements. 



Detector 


Performance 


Coverage 


Vertex detector 


a b = 7.0 © (20.0/p) / sin 3/2 9 /mi 




cos# 


< 0.90 


Central drift chamber 


a PT /pT = 1-1 x 10~V © 0.1 % 




cos 9 


< 0.95 


EM calorimeter 


a E /E = 15 % /y/E © 1 % 




cos 9 


< 0.90 


Hadron calorimeter 


a E /E = 40 % /VE © 2 % 




cos 9 


< 0.90 



Table 1: ACFA study parameters of the JLC detector, where p, pt, and E are measured in 
units of GeV. 



3 Event Selection 

3.1 Basic Reconstruction Strategy 

Since the top quark decays almost 100% into a b quark and a W boson, the signature of 
a ti production is two b jets and two W bosons in the final state. These W bosons decay 
subsequently either leptonically into a lepton plus a neutrino or hadronically into two jets. 
According to how the W bosons decay, therefore, there will be three modes of final states: (1) 
six jets, where both of the W's decay hadronically, (2) one lepton plus four jets, where one 
of the W's decays leptonically and the other hadronically, and (3) two leptons plus two jets, 
where both of the W's decay leptonically. 

In order to reconstruct the momentum vector of the top quark, we will use the lepton- 
plus-4-jet mode, for which we can reconstruct the t(i)-quark momentum as the momentum sum 
of the b(b) jet and the two jets from the hadronically-decayed W + (W~), while we can tell the 
charge of the hadronically-decayed W from the charge of the lepton. 

In the lepton-plus-4-jet mode, two of the four jets are b(b) jets directly from the t(t) 
quarks, while the other two are from the W boson that decayed hadronically. Therefore, if one 
can identify the b and b jets, remaining two jets can be uniquely assigned as decay products 
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of the W boson. The other W boson can be reconstructed from the lepton and the neutrino 
indirectly detected as a missing momentum. Remaining task is then to decide which b(b) jet 
to attach to which VT-boson candidate, in order to form t(t) quarks. Since the t(t) quarks 
are virtually at rest near the threshold, a b(b) jet and the corresponding W boson fly in the 
opposite directions. We can thus choose the correct combination by requiring the b(b) jet and 
the W boson be approximately back-to-back. 

In reality, however, 6(6)-quark tagging is not perfect and can be performed only with 
some finite efficiency and purity: there could be more than two &(6)-jet candidates in a single 
event. In addition, b and b quarks can be emitted in the same direction. In such a case, a 
wrong combination could accidentally satisfy the back-to-back condition. These facts sometimes 
prevent us from uniquely assigning each jet to its corresponding parton, resulting in multiple 
solutions for a single event. Moreover, the leptonically-decayed W is poorly reconstructed in 
practice, since the neutrino momentum is strongly affected by ISR, beamstrahlung, as well as 
other possible neutrinos emitted from the b(b) jets. In order to overcome these difficulties, we 
will need some sophisticated method. We defer discussion of such a method to the next section 
and examine here the extent to which the aforementioned basic reconstruction strategy works. 



3.2 Event Selection Procedure 

The lepton-plus-4-jet-mode selection started with demanding an energetic isolated lepton: Eg > 
18 GeV and £a4° C one < 18 GeV, where Eg is the lepton's energy and -Ei4° COI1 e is the energy sum 
of particles within a cone with a half angle of 14° around the lepton direction excluding the 
lepton itself.^] When such a lepton was found, the rest of the final-state particles was forced 



clustering to four jets, using the Durham clustering algorithm p3| . Two-jet invariant mass was 
then calculated for each of the six possible combinations and checked if it was between 65 GeV 
and 95 GeV, in order to select a jet pair which was consistent with that coming from a VT-boson 
decay. For such a jet pair the remaining two jets, at the same time, had to be identified as b(b) 
jets, using flavor tagging based on the impact parameter method. The hatched histogram in 
Fig. H| is the 2-jet invariant mass distribution of all the possible pairs out of the four jets, while 
the solid histogram being that with the 6-tagging. It is seen that this procedure dramatically 
improved the purity of the W boson sample. It should also be stressed that these selection 
criteria are very effective to suppress background processes such as e + e~ — > W + W~ and provide 
us with an essentially background-free ti event sample. 

The remaining task is to decide which b(b) jet to associate with which W candidate. For 
a b(b)-jet candidate, the right W boson partner was selected by requiring the back-to-back 
condition as described above. Fig. is a scatter plot of the acoplanarity angles of the two 
possible b-W systems where horizontal and vertical axes are the angles of b-Wg u and b-W2-j e t 
system, respectively. b-W pairs having 9 acop ^-w) < 60° was regarded as daughters of the t(t) 
quarks. 

The selection efficiency after all of these cuts was found to be 15% including the branching 
fraction to the lepton-plus-4-jet mode of 29%. 



4 Kinematical Fit 

The event selection described above yields a very clean ti sample. As noted above, however, 
the sample is still subject to combinatorial backgrounds, if we are to fully reconstruct the final 

* The Ei cut was chosen to be the kinematical limit for the lepton from the W — > Iv decay. On the other 
hand, the cone-energy cut was optimized to achieve high purity, while keeping reasonable efficiency. 
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Figure 1: Invariant mass distributions of the 2-jet systems reconstructed as VT-boson candi- 
dates. Hatched and solid histograms correspond to before and after double ft-tagging, respec- 
tively. The locations of the W mass cuts are indicated with arrows. 
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Figure 2: Scatter plot of the acoplanarity angles corresponding to two b-W systems, where 
horizontal and vertical axes are angles of b-Wi v and b-W 2 -j e t systems, respectively. 
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state by assigning each jet to a corresponding decay daughter of the t or t quark. We thus need 
a well-defined criterion to select the best from possible multiple solutions. It is also desirable 
to improve the measurement accuracies of those kinematical variables which are suffering from 
effects of missing neutrinos (such variables include momenta of b, b or the neutrino from a W 
itself). 

The ti system produced via e + e~ annihilation is a heavily constrained system: there are 
many mass constraints in addition to the usual 4-momentum conservation. At e + e~ linear 
colliders, thanks to their well-defined initial state and the clean environment, we can make full 
use of these constraints and perform a kinematical fit to select the best solution and to improve 
the measurement accuracies of the kinematical variables of the final-state partons. 

4.1 Parameters, Constraints, and Likelihood Function 

For the lepton-plus-4-jet final state, there are 10 unknown parameters to be determined by the 
fit: the energies of four jets, the 4-momentum of the neutrino from the leptonically-decayed W 
boson, and the energies of the initial-state electron and positron, provided that the jet directions 
as output from the jet finder are accurate enough, the error in the 4-momentum measurement 
of the lepton from the leptonically-decayed W can be ignored, and that the transverse momenta 
of the initial-state electron and positron after beamstrahlung or initial-state radiation or both 
are either negligible or known from a low angle e + /e~ detector system^] (see Fig. ^). 




Figure 3: Schematic diagram showing parameters and constraints relevant to the kinematical 
fit described in the text. The boxed parameters are unknown and to be determined by the fit. 

The requirements of 4-momentum conservation and the massless constraint for the neu- 
trino from the leptonically-decayed W reduce the number of free parameters to 5. We choose, 
as these free parameters, the energies of the four jets and the initial longitudinal momentum 
(the difference of the energies of the initial-state electron and positron). 

These five unknown parameters can be determined by maximizing the following likelihood 
function: 

L = (6 %( E r n > E f)) ■ P r w+ ■ Pr w - ■ Pr« ■ ^/I, (1) 
/=1 

where Pg is a resolution function for jet / and is Gaussian for / = 1 and 2 (jets from the 
hadronically-decayed W) as given by the detector energy resolution. For / = 3 and 4 (jets 

t In addition, there will be some finite transverse momenta due to a finite crossing angle of the two beams. 
These transverse momenta are, however, known and can be easily incorporated into the fit. 



5 



from the b and b quarks) the resolution function is the same Gaussian convoluted with the 
missing energy spectrum due to possible neutrino emissions. For the two W bosons in the final 
state, we use a Breit-Wigner function P^ w instead of 5-function-like mass constraints. is a 
weight function coming from ISR and beamstrahlung effects. This distribution was calculated 
as a differential cross section as a function of the energies of initial-state electron and positron, 
taking into account the tt threshold correction as described in Sec. ||. 

The remaining factor, Pr t{ , controls the mass distribution of the t and i quarks and 
has been introduced to take into account the kinematical constraint that the t and i cannot 
be simultaneously on-shell below threshold (see Fig. f| which shows Pp tJ distribution below 
tt threshold). Pr ti distribution is a dynamics- independent factor which is extracted from the 
theoretical formula for the threshold cross section. 
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Figure 4: P rtf distribution below tt threshold. 



4.2 Results 

We performed the maximum likelihood fit for the selected sample. The maximum likelihood 
fit provided us with a well-defined clear-cut criterion to select the best solution, when there 
were multiple possible solutions for a single event: we should select the one with the highest 
likelihood. 

Figs. |5|-a) and -b) are the reconstructed W mass distributions for the leptonically and 
hadronically-decayed W bosons, respectively, before (hatched) and after (solid) the kinematical 
fit. The figures demonstrate that the Breit-Wigner factors (Pr w ±) m the likelihood function 
properly constrain the W masses as intended. 

Fig. ||-a) plots the reconstructed mass for the t(t) decayed into 3 jets against that of the 
i(t) decayed into a lepton plus a b jet, before the kinematical fit. The strong negative correlation 
is due to the fact that the neutrino from the leptonically-decayed W is reconstructed as the total 
missing momentum. Figs. |^-b) and -c) are the projections of Fig. ^|-a) to the horizontal and 
vertical axes, respectively, showing systematic shifts of the peak positions.^ Figs. |§-d) through - 
f ) are similar plots to Figs. |B|-a) through -c) after the kinematical fitting, while Figs, ^-g) through 
-i) are corresponding distributions of generated values (Monte-Carlo truth). The kinematical 
fit sent most of the events to the L-shaped region indicated in Fig. P-d), as it should, and made 

$ This is in contrast with the result in Q , where a quite tight set of cuts was imposed upon the reconstructed 
W and t masses, and consequently their peak shifts were less apparent at the cost of significant loss of usable 
events. The goal of this study is to establish an analysis procedure to restore those events which would have 
been lost, by relaxing the tight cuts while keeping reasonable accuracy for event reconstruction. 
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Figure 5: Reconstructed W mass distributions for (a) leptonically and (b) hadronically-decayed 
W bosons, before (hatched) and after (solid) the kinematical fitting. 

the distribution look like the generated distribution shown in Fig. |6|-g). Consequently, the peak 
shifts observed in the Figs. ^|-b) and -c) have been corrected as seen in Figs. ||-e) and -f). There 
are, however, still some small fraction of events left along the minus 45° line. These events 
were so poorly measured that it was impossible to restore. The cut (angled region) indicated in 
Fig. |6|-d) allowed us to remove them without introducing any strong bias on the reconstruction 
of the kinematical variables. 

Now the question is how the above constraints improve the parameters of the fit such as 
the energies of b and b jets, the direction and the magnitude of the missing neutrino from the 
leptonically-decayed W, on which we expect significant influences. Figs. |7|-a) and -b) plot the 
difference between the reconstructed and the generated energies of the b (b) quark attached 
to the leptonically-decayed W and that of the b (b) attached to the hadronically-decayed W, 
respectively, before (hatched) and after (solid) the kinematical fit. The plots demonstrate that 
the kinematical constraints recover the energies carried away by neutrinos from the b or b decays. 
The improvement is more dramatic for the direct neutrino from the leptonically-decayed W, 
which is reconstructed as the total missing momentum; see Figs. 0-c) and -d) which show 
distributions of the difference of the reconstructed and generated neutrino energies (AE V ) and 
directions (A9 V ). 

The improvements in these kinematical variables are reflected to the improvements in the 
reconstructed W energies and directions as shown in Fig. §-a) for the energy of the leptonically- 
decayed W, -b) for the hadronically-decayed W, and -c) for the direction of the leptonically 
decayed W. We can see dramatic improvements in all of these distributions, although the 
improvement in the direction of the hadronically-decayed W was less dramatic. 

Finally, we will examine the effects of the kinematical fit on the measurements of the 
direction and the magnitude of the top quark momentum. In Figs. ^|-a) and -b), the difference 
of the reconstructed and generated directions of the t or i quark is plotted against the generated 
top momentum, before and after the kinematical fit, respectively. We can see appreciable 
improvement by the fit. Nevertheless, since the top quark direction becomes more and more 
difficult to measure as the top quark momentum decreases, the resolution is still somewhat poor 
in the low momentum region. The angular resolution is largely determined by the reconstruction 
of the t or i decayed into 3 jets. Remember that the resolution improvements were less significant 
for the hadronically-decayed W, since the power of the constraints was used up mostly to recover 
the momentum information of the direct neutrino from the leptonically-decayed W and the 
energy resolution for jets from the W was left essentially unimproved. The improvement in the 
measurements of the top quark direction is mostly coming from the improvement in the b or b 
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Figure 6: Scatter plot of the reconstructed t(t) mass from 3 jets versus that from a lepton plus 
a fo-jet (a) before the kinematical fit, together with (b) its horizontal/5£z/ and (c) vertical/3-jet 
projections, (d) through (f) are similar plots after the kinematical fit and (g) though (i) are 
corresponding plots for generated values. 



S 





Figure 7: Distributions of the difference of the reconstructed and generated energies of the b or 
b jet attached to (a) leptonically-decayed and (b) hadronically-decayed W bosons, together with 
distribution of the difference of the reconstructed and generated (c) energies and (d) directions 
of the direct neutrino from the leptonically-decayed W , before (hatched) and after (solid) the 
kinematical fit. 




Figure 8: Distributions of the difference of the reconstructed and generated energies of (a) 
leptonically-decayed and (b) hadronically-decayed W bosons, and (c) distribution of the dif- 
ference of the reconstructed and generated directions of the leptonically-decayed W, before 
(hatched) and after (solid) the kinematical fit. 
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jet measurement. By the same token, the effect of the fit on the measurement of the magnitude 
of the top quark momentum is also less dramatic compared to that on the leptonically-decayed 
W. 
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Figure 9: The difference of the reconstructed and generated directions of the t or t quark 
plotted against the generated top momentum, (a) before and (b) after the kinematical fit. 

In the case of the 6-jet mode, for which there is no direct energetic neutrino from W's, we 
can use the power of the constraints to improve the jet energy measurements. Consequently, 
we may expect more significant improvement in the top quark momentum measurement. 



5 A Possible Application 

We discuss a possible application of our kinematical reconstruction method. Let us consider 
measurements of the decay form factors of the top quark in the ti threshold region. We assume 
that deviations of the top-decay form factors from the tree-level SM values are small and 
consider the deviations only up to the first order, i.e. we neglect the terms quadratic in the 
anomalous form factors. Then the cross sections depend only on two form factors /f and 
in the limit m;, — > although the most general tbW coupling includes six independent form 



factors |26| : 



u f Lp ^ u Pw, fRp • 



u{p t ), (2) 



where P L = (1 - 7 5 )/2 and P R = (1 + 7 5 )/2. At tree level of the SM, f[ = 1 and = 0. 
A variation of changes only the normalization of the differential decay width of the top 
quark, whereas a variation of changes both the normalization and the shape of the decay 
distributions. Thus, we expect that the kinematical reconstruction is useful for disentanglement 
of the two form factors and in particular for the measurement of f^. For simplicity we assume 
fi=l hereafter .| Since transverse W (denoted as Wt) is more sensitive to /vf than longitudinal 
W (Wl), our strategy is to extract Wt using the angular distribution of W (in the rest frame 
of t) and the angular distribution of £ (in the rest frame of W). It is well known that Wt 
is enhanced in the backward region cosOw — —1, where the angle 9w is measured from the 
direction of the top quark spin in the t rest frame. Also, we may enhance Wt by collecting 
£ emitted in the backward direction cos#£ ~ —1, where the angle 8g is measured from the 
direction of — p t in the W rest frame. These features are demonstrated in Figs. 0: We plot[] 

§ In order to determine simultaneously, we may, for instance, use independent information from the 
measurement of the top width |tJ . 

^ We used the helicity amplitudes given in [E6| for calculating these differential decay widths. 
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(a) the differential decay width for the decay of the top quark with a definite spin orientation 
dT(ti — > htv) I \d cos 9\yd cos 9g) for f^ = and (b) the difference of the differential widths for 
f2 = 0.1 and for = 0. The plots show that we may measure f^, for instance, from the 
ratio of the numbers of events in the regions cos^vk, cos9e < and cos 9w, cos 9e > 0. Since 
a most prominent result of our kinematical fitting is the improvement in the measurement of 
the momentum of leptonically decaying W, we expect that the sensitivity to /<f would increase 
after the kinematical fitting.^ 




Figure 10: (a) Normalized differential decay width Af~ 1 dr(t-^ — > htv) j{d cos 9y/d cos 9 'i) for 
f% = 0- O 3 ) Difference of the normalized differential decay widths for /g = 0.1 and for f% = 0. 
In both figures the differential widths are normalized by M = T t x Br(W / — > tv) for f% = 0. 



It is quite advantageous to investigate decay properties of the top quark in the tt threshold 
region as compared to the open-top region E ^> 2m t because of several reasons. It is useful 
that the top quark can be polarized close to 100% in the threshold region j27], |28j. This is 



clear in the above example. Furthermore, we are almost in the rest frame of the top quark. 
In the above example, the top quark is highly polarized in its rest frame. Hence, the event 
rate expressed in terms of cos 0w and cos 6^ is a direct measure of the amplitude-squared, 
| J2i=L,T -A(ti — > hWi) x A(Wi — > tv)\ 2 (without phase-space Jacobian), which allows for simple 
physical interpretations of event shapes. On the other hand, we do not gain resolving power for 
the decay form factors by raising the cm. energy. This is in contrast with the measurements 
of the tt production form factors. 



6 Summary and Conclusions 

To make maximum use of future e + e~ linear colliders' experimental potential, the top quark 
reconstruction in the lepton-plus-4-jet mode has been studied under realistic experimental con- 
ditions of e + e~ — > tt process near its threshold. As a new technique to fully reconstruct ti 
final states, we have developed a kinematical fitting algorithm which aims to reconstruct the 
momentum- vector of top quark more accurately. 

The missing energy carried away by neutrinos from bottom quark decays has been recov- 
ered by the kinematical fitting. However, the effects of the kinematical fitting on the top quark 
momentum are not as dramatic as we wanted. This is because the top quarks are almost at rest 

II We can measure f% also from the energy distribution of I. (The angular distribution of I is insensitive to 
f% [HI Hi-) ln this case, we do not need to reconstruct the W momentum. From a rough estimate, however, 
we find that the method using the W and I angular distributions has a higher sensitivity to f^. 
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in the threshold region and therefore their momenta are difficult to measure. Moreover, in the 
lepton-plus-4-jet mode many constraints are used up by recovering the information on the neu- 
trino from leptonically-decayed W bosons. On the other hand, the remarkable improvements 
of the energy resolution of fe-jets and the angular and energy resolutions of leptonically-decayed 
W's have been achieved by the kinematical fitting. These improvements should benefit the 
form factor measurements in general. As a possible application, we considered measurements 
of decay form factors including on which correct reconstruction of the leptonically-decayed 
W may have a large impact. 

There have been a number of theoretical studies on measurements of the top-quark pro- 
duction and decay form factors using the e + e~ — > ti process. Many of these analyses assumed 
either the most optimistic case or the most conservative case with respect to the kinematical 
reconstruction of event profiles. In the former case, one assumes that the momenta of all the 
particles (including t and W) can be determined precisely; in the latter case, one uses only par- 
tial kinematical information, e.g. the direction of b, and the energy and momentum of t. Our 
analysis indicates that both assumptions are not the most realistic approximations in actual 
experimental situations. 

As further studies, we will move on to investigate expected sensitivities to various anoma- 
lous couplings in the lepton-plus-4-jets mode systematicaly. We also plan to extend the kine- 
matical fitting algorithm to the 6-jet mode, where we expect more significant improvements of 
the momentum resolution of the top quark and combinatorial background. 
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